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Abstract 

The glob model of Lichard and Van Hove and the modified soft annihilation 
model (MSAM) of Lichard and Thompson are used as a phenomenological 
tool for relating results from various experiments on soft photon production in 
high energy collisions. The total phenomenological expectation is composed of 
contributions from classical bremsstrahlung, the soft annihilation model and 
the glob model. The empirical excess above the background from hadronic 
decays at very small longitudinal momenta of photons is well reproduced, 
as well as that for transverse momenta pt ^ 10 MeV/c. Some data do 
not require the glob model and MSAM components in the phenomenological 
mixture, but do not exclude them. On the basis of consistency of all data with 
the total theoretical expectation we argue that the results of all experiments 
are mutually consistent. The models are unable to describe the excess of 
ultrasoft photons (pr 10 MeV/c), seen by some, but not all, experiments. 
This may indicate an as yet unknown projectile-mass-dependent production 
mechanism. Possible relations of soft photon production to other phenomena 
are discussed. A simple-to-use, but physically equivalent version of the glob 
model is developed, which enables an easy check of presented results. 
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I. INTRODUCTION 



The inclusive production of soft photons in particle collisions at high energies has been 
studied in several experiments [[L]-[| using various experimental techniques. 

In a pioneering bubble chamber experiment performed at the Stanford Linear Accelerator 
Center (SLAC), Goshaw et al. measured the transverse momentum distribution of photons 
produced in ir + p collisions at 10.5 GeV/c. Using their own data about hadronic final states 
they showed that the bremsstrahlung formula of classical electrodynamics |I(J was able to 
account for all the observed photon yield. 

A few years later, the experimental group led by Yves Goldschmidt-Clermont inves- 
tigated the photon production in 70 GeV/c K + p collisions using the Big European Bub- 
ble Chamber (BEBC) at the CERN Super Proton Synchrotron (SPS). In contrast to the 
observation of the SLAC experiment they found many more photons with very small lon- 
gitudinal momentum {2\p L \/^fs < 0.005) than they expected on the basis of the classical 
bremsstrahlung formula and their charged hadron production data. Similarly, the unpub- 
lished fill] photon transverse momentum spectrum from the same experiment exhibited a 
significant excess over the bremsstrahlung estimate for pt 60 MeV/c. 

The results obtained by the Axial Field Spectrometer (AFS) collaboration at the CERN 
Intersecting Storage Rings (ISR) H were in conformity with both previous experiments. 
They neither required an additional source to explain the observed photon signal, nor ex- 
cluded an anomalous component of the size indicated by the BEBC experiment. But they 
did rule out a strong increase of anomalous photon production with rising collision energy 
(the invariant energy a/s in the AFS experiment was 63 GeV, compared to 11.5 GeV in the 
BEBC experiment). 

In a subsequent AFS experiment J|], the direct photon production was studied in pp and 
aa collisions at ^ = 63 GeV. The photon momentum range explored (0.1 GeV/c < pt < 
1 GeV/c) does not overlap the range in this paper (j>t ^0.1 GeV/c). This is the reason 
why we will not include these data in our comparison. Let us only note that the results of 
Q did not show any excess beyond the limits dominated by systematic uncertainties. 

The inclusive yield of photons from deep inelastic \ip scattering at 200 GeV was measured 
by the European Muon Collaboration (EMC) at the CERN SPS ||. After subtracting the 
contributions from hadron electromagnetic decays and Bethe-Heitler muon bremsstrahlung, 
residual photons were observed at a mean level of 0.15 ± 0.06 per event. 

The first results of the HELIOS (NA34) experiment at the CERN SPS [|] on soft photon 
production in p-Be and p-Al interactions at 450 GeV/c indicated a marked excess in the 
region 4 MeV/c < px < 20 MeV/c. In addition, the yield seemed to increase approximately 
as the square of the associated hadron multiplicity, which would have signified a collective 



production mechanism [[T2]| . These findings were not confirmed by a more extensive later 
experiment || by some members of the same group. 

In the experiment performed by the EHS-NA22 Collaboration |7j at the CERN SPS, 
the European Hybrid Spectrometer (EHS) was equipped with the Rapid Cycling Bubble 
Chamber (RCBC) as vertex detector. The inclusive soft photon production was studied in 
n + p and K + p collisions at 250 GeV/c. The results confirmed the existence of an anomalous 
prompt photon signal with very similar properties as seen in the BEBC experiment |p|Jil~||. 
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The SOPHIE/WA83 experiment ||, which was performed at the CERN SPS, used the 
OMEGA spectrometer supplemented with two electromagnetic calorimeters. It differed in 
three important respects from the previous two hadronic collision experiments that saw 
the anomalous soft photon signal. First, it used a beam of negative particles (pions with 
momentum of 280 GeV/c). Second, it was a purely electronic experiment without utilizing 
a bubble chamber. Third, it explored a kinematic region where the contribution of gammas 
from hadronic decays is relatively small. The fact that the results of this experiment, as 
concerns the existence and approximate magnitude of the anomaly, are in conformity with 
those of previous ones ||[7J is therefore especially valuable. 

Antos et al. || used a modified setup of the HELIOS experiment || and measured the in- 
clusive pt spectra of soft photons produced at central and slightly backward rapidities in 450 
GeV/c p-Be collisions. Two independent photon measurement methods with corresponding 
detectors and analysis chains were used in parallel: (i) a combination of gas chambers, con- 
verter plates and a BGO matrix for the conversion method; (ii) a BaF 2 array for the time 
of flight photon identification method. The authors of observed a significant excess of 
direct photons at very low pr (< 15 MeV/c) above the background from hadronic decays. 
This excess was consistent with the expected contribution from hadronic bremsstrahlung, 
calculated from the classical electrodynamics formula |TTJ . 

On the theoretical side, the results of the EMC photon experiment H were described 



|13| using the soft annihilation model ||i4|,|15| . Some discrepancy between experimental and 
theoretical spectra stimulated the creation of the modified soft annihilation model (T^|, which 
allowed, contrary to the original model, also gluons in the intermediate parton state. We 
will present this model in some detail in Sec. [TIT]. Every attempt to explain the photon 
excess seen in the BEBC experiment 0] on the same footing remained unsuccessful. 

To our knowledge, the first theoretical article that addressed the issue of very soft anoma- 
lous photons || was that by Van Hove |I7| . He offered a common explanation of several 
ultrasoft phenomena observed experimentally in central rapidity regions of high energy 
collisions. These phenomena are characterized by very low transverse momenta (produc- 
tion of soft pions or photons) or very short distances in rapidity (fluctuations in rapidity 
distribution- "intermittency"). Van Hove argued that they could be regarded as manifesta- 
tions of the occurrence, in at least some of the collisions, of an intermediate parton system 
with considerable lifetime and spatial extension. This system carries only a part of the 
collision energy and momentum and the rest of the event is "standard" . The most natural 
framework for understanding the production mechanism and properties of such a system 



is provided by the QCD parton shower model [|T8|. It is usually assumed that the shower 
development stops when parton virtualities fall to Qq ~ 1 GeV. The partons with this vir- 
tually are supposed to enter the hadronization process. Van Hove proposed that some of 
them may continue in showering, producing a large multiplicity (N p ~ 30) system of very 
soft partons, a glob of cold quark-gluon plasma. The intermediate parton system consists 
of one or several globs with masses Mq w 1 GeV/c 2 . As discussed in [[17]], globs need much 



time to hadronize, because hadronization requires a very drastic rearrangement of partonic 
wave functions. During their long lifetime, globs can produce photons in the subprocesses 
q + q — > 7 + g and g + q(q) — > 7 + q(q). A quantitative model for soft photon production 
based on the Van Hove ideas |L7j was constructed and successfully compared to the BEBC 
data P,|TT| by him and the present author in . 
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Barshay investigated the role of pion condensation f21"|j22}1 in soft photon production. 
He showed that this collective coherent mechanism implied that the photon emission was 
proportional to the square of the associated pion multiplicity. At that time, this feature 
seemed to be indicated by the preliminary HELIOS data ||. 

Shuryak |23| showed that the backward reflection of pions at the boundary of a hadronic 
system, induced by the modification of the pion dispersion curve by many-particle interac- 
tions, strongly increases soft photon emission. A complex view of the dense interacting pion 
matter was presented in the subsequent work p4| together with the implications for exper- 
imentally observed phenomena. The photon yield was calculated by applying the classical 
bremsstrahlung formula along the paths of many times rescattered pions. As stated in the 



original paper |24]| , this approach could not describe the observed excess of real photons at 
low transverse momenta. 

V. Balek, N. Pisutova, and J. Pisiit |25| (a more detailed discussion of some issues can be 
found in their later papers: in collaboration with Zinovjev, and P7|1 ), first briefly reviewed 
experimental information on, and theoretical understanding of the production of very soft 
photons. Then they studied some effects which were not considered in the bremsstrahlung 
calculation of the HELIOS Collaboration || and found them very small. They also suggested 
two photon production mechanisms: 1) shock waves in the system of final state hadrons, 
2) bremsstrahlung emitted by a quark that tries to escape from the intermediate parton 
system. They concluded that none of them could explain the soft photon anomaly observed 
in experiments IfJH . 

Czyz and Florkowski [28] used the classical bremsstrahlung formula to calculate soft 



photon emission within the framework of the boost-invariant color-flux tube model p9 



The ratio of their photon production rate to the classical one depends strongly on the 
direction of the photon and is very sensitive to the assumed mass of quarks. For photons 
emitted perpendicular to the collision axis and m q ~ 10 MeV/c 2 it can reach 10. The strong 
angular dependence is a feature that saliently distinguishes this model from the essentially 



isotropic glob model [19 



Only a few of the theoretical approaches mentioned above have attempted to make 
a detailed comparison with actual experimental data, including the very important and 
delicate matter of experimental acceptance and cuts. 

While technical and experimental problems are substantial, an important problem in 
assessing the reality of the anomalous photon production is that the experiments have been 
performed under very different conditions. The varied experimental conditions have included 
different projectiles and targets, different collision energies, different instrumental setups 
with different photon momentum coverage and acceptances. So it is very difficult to say 
whether they are consistent among themselves or not, whether they witness about the same 
phenomenon, whether their contradictory claims are really significantly inconsistent. In this 
complicated situation we see only one way of pursuing the matter: to use some theoretical 
model as a tool for connecting various experiments. The guiding idea is that if a model is 
able, after modulated by the experimental acceptances, to describe several pieces of data 
without tuning its parameters for each particular case, then the data are mutually consistent. 
It may also indicate that the model is physically sound, but this need not always be the 
case. The model may simply simulate an important conventional contribution which was 
not properly taken into account in every experimental analysis considered. 
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The purpose of the present paper is to assess the consistency of various data on soft 



photon production in hadronic interactions by using the glob model [19] and the modified 



soft annihilation model [I6[] as interpolation tools. 

In the next section, we deal with the glob model of Van Hove and the present author. 
We first recapitulate its assumptions and equations, as well as the way its parameters were 
chosen. Then we develop a simpler, but physically equivalent mutation of the glob model, 
which will be used throughout this paper. It can easily be implemented by anyone wishing 
to perform her or his own calculations. Section |T| deals similarly with the modified soft 
annihilation model of Thompson and the present author. Unfortunately, here we are unable 
to offer a user-friendly version. The central part of the paper is Sec. [TV], where we show the 
results of model calculations and their comparison with experimental data. Conclusions are 
summarized and commented upon in Sec. M, where we also discuss the possible relation of 
anomalous soft photon production to other phenomena. 



II. GLOB MODEL 

A. The original version 

As mentioned in Sec. |, we assume that in some percentage of high-energy inelastic 
collisions, a long-living, large, and dense system consisting of light quarks, antiquarks, and 
gluons is formed. For the physical parameters of such a system, called glob, we will use the 
original estimates |T7|P^[ : the glob mass Mq = 1 GeV/c 2 ; the number of partons within 
a glob N p = 40; the gluon/parton number ratio = 0.5; and the light quark number ratio 
u/d = 1. 

Because of the fixed invariant energy, the momentum distribution of partons is governed 
by the microcanonical distribution. The mean number of photons emitted by a glob is 
therefore given by the expression 



i Np 



1/2 



x 'yv^j 



-S\J2E k -M\5 £ Pfe Il-ir 1 . (2- 1 ) 

\k=l J \k=l J fc=l E k 

where is the momentum of the kth parton in the glob rest frame and &ij(sij) is the cross 
section for photon production in head-on collisions of the zth and j'th parton. The processes 
q + q — > 7 + g and g + q(q) — > 7 + q(q) are considered in the lowest order of QED and QCD. 
The phase-space normalization constant is 

(N p \ / N p \ N p ,3 

EvmU Ep* iIt- (2 - 2) 
k=l I \fc=l / k=l & 



The total cross section cr^ in Eq. (|2.1|) can be written as 

dipij daij 



Vij(sij) = I dUj ^ df J , (2.3) 



Kl 
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where tij is the four-momentum transfer squared from one parton to the photon in the 
collision of the ith and jth partons and ifiij is the azimuthal angle of the photon in their 
rest frame. After inserting ( |2.3| ) into ( |2.1| ) and mapping the independent variables onto a 
(3N p — 2)-dimensional unit cube ||30|| , we arrive at the equation that served in [1J| as the 



master equation for a Monte Carlo generator of the photon momenta in the glob rest frame. 

To transform the photon momentum from the glob rest frame to the collision center-of- 
mass frame, we need to know the glob momentum distribution in the latter. In (jTJJ it was 
assumed that the distribution in the glob rapidity and transverse momentum squared 

cPN 

, A = N G F G (y G , Pt,g) > (2-4) 
dy G dp 2 TG V - J 

with the function F G normalized to unity, factorizes. To proceed further, a gaussian glob ra- 
pidity distribution with the width proportional to the maximum cms rapidity (i.e. logs) was 
chosen, with {y G ) l l 2 = 0.6 for ^fs = 11.51 GeV (the BEBC experiment 0). The exponential 
p\ G distribution was assumed to be collision energy independent with (pt,g) — 0.3 GeV/c. 
In each Monte Carlo "event", the rapidity and transverse momentum were generated, in- 
creasing the number of independent variables to 3N P . Then the photon momentum was 
transformed into the collision center-of-mass frame. 

To complete the model, the mean number of globs per event N G is multiplied by the 
mean number of photons per glob to give the mean number of photons per event n 7 . 
Keeping in mind that the cross sections of photon production subprocesses are proportional 
to the strong coupling constant a s , the overall multiplication constant that fixes the absolute 
normalization of photon yield is 

B = a s N G ^r • (2.5) 

In order to fit the BEBC data g], its value was fixed in |B| at B = 5 x 10~ 3 c -1 fm~ 2 . 
To get the predictions for photon spectra at different collision energies or with different 
incident particles it was suggested in |1| to scale the mean number of globs N G according 



to empirically known hadron multiplicity. 

The first numerical realization of the glob model [19], which we have just briefly sketched, 
had serious disadvantages. With N p fixed at 40, the dimension of the integral which was to be 
evaluated by the Monte Carlo method was 120. Determination of a detailed distribution in 
photon momentum with low statistical errors thus became a computer-time-consuming task. 
Also the computer code was complicated and difficult to use. It was practically impossible to 
use the glob model as an event generator for Monte Carlo studies of experimental setups. In 
the next subsection we remove these drawbacks of the original glob model |19|] and suggest 
a computing scheme which is easy to reproduce. 



B. A simplified version of the glob model 

In the glob model it is assumed that the internal properties of the glob (total invariant 
energy of partons M G c 2 , their number N p , as well as their momentum and flavor distribu- 
tions) depend neither on the type of process under study nor on the incident energy. This 
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allows us to construct a new version of the glob model, which is physically equivalent to the 
original one, but much easier to use. 

The invariant photon distribution in a general reference frame can be written as a con- 
volution of the photon distribution in the glob rest frame with the momentum distribution 
of globs in the general frame 



(i 3 n 7 f 3 d 3 N G 



n 



a 



7 <i 3 p 7 J d PG ^3p G 7 ^3p* ( 6) 



The asterisk refers to quantities in the glob rest frame and 



According to the basic assumptions of the glob model, the photon momenta are distributed 
isotropically in the glob rest frame. It is thus sufficient to consider only the energy spectrum 
in the latter frame. We can write 

1 ~ • (2-8) 



A study of this quantity in the framework of the original Monte Carlo version of the glob 
model has showed that the function f(E*) can be parametrized as 

f(E;)=aexp{-bE;-c(E;) 2 } , (2.9) 

with parameters a, b, and c given in Table |. Comparison of the photon energy distribution 
in the glob rest frame based on the parametrization ( |2.9p with the original Monte Carlo 
calculation in Fig. [I] shows that they are identical. 

Introducing the distribution in the glob rapidity and transverse momentum squared 
and using Eq. (|2.8|) , we can cast Eq. ( |2.6| ) in the form 

<i 2 n 7 B p Tn 



dPL,ydp T ,~f 4:71 E y 

where B is again given by Eq. (|2.5|) and 



J E G (y G , p 2 TjG ) f (e;) dy G dp^ G d(j) , (2. 10) 



E* = —— (E G Ej - p LtG p Ltl - Pt,gPt,~< cos0) . (2.11) 
M G 



Even if the integrals in Eq. Q2.10] ) can be evaluated by more conventional numerical methods, 
a Monte Carlo approach is convenient to utilize, especially if the experimental acceptance 
has to be taken into account. 



III. MODIFIED SOFT ANNIHILATION MODEL 



The original soft annihilation model (SAM) [[14] was inspired by Bjorken and Weisberg 



[511 who suggested an explanation of an unexpectedly large production of low- mass lepton 
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pairs discovered in many experiments f32|| . In SAM, dileptons arise from an intermediate par- 
ton system (IPS) by annihilation of quarks and antiquarks (qq — > and qq — > gluori) 
produced in the initial stages of the reaction. Unlike the IPS system of Van Hove [117]], the 
IPS of SAM carries all the collision energy and is the only source of final state hadrons. Its 
parameters are therefore fixed by data on the production of hadrons |[33|| . In 1981, SAM 
was compared |L5| to all the dilepton data available at that time and appeared to be in 
satisfactory agreement with them. 

Later on, the data on electron production at the CERN Intersecting Storage Rings (ISR) 
were extended to lower values of the transverse momentum [ 31 ] . The SAM was not able to 
follow a steep rise of the e + /it ratio with decreasing transverse momentum. The comparison 
of SAM (with the subprocess qq — > 7 + g) to the EMC photon data was not without flaws 
either. Both longitudinal and transverse spectra (in the hadronic system rest frame) were 
flatter than the experimental one [13|]. To bring the calculations closer to reality, a gluon 
component was introduced into the IPS in the modified soft annihilation model [|16| together 
with corresponding subprocesses g + q(q) — > 7 + q(q) and g + q(q) — > + q(q). With 
the quark and antiquark multiplicities fixed at the same values as in SAM (in order not to 
change outgoing hadron multiplicities), the mean energies of partons became smaller, which 
also made the spectra of dileptons and photons softer. This improved the agreement not 
only with the results of and pijl , but also with some older dilepton data (see flTB]] ). We 
refer the reader to [16] for technical details of the model, which we will use also in this work. 
Let us only note that the model, when applied to soft photon production, does not contain 
any free parameters. For each projectile, target, and incident energy combination, the mean 
number of partons in IPS was fixed by the mean hadron multiplicity. 



IV. COMPARISON OF MODELS TO DATA 

In this section we compare the outcomes of the glob model Jl9| and the modified soft 



annihilation model (MSAM) [0 with existing data. For data which provide absolute nor- 
malization, the procedure is straightforward. They are those listed in the Introduction 
except for the AFS and Antos et al. [0] experiments. 

In paper || , the AFS Collaboration presented the pT-spectrum of photons observed in 
pp collisions at \/s = 63 GeV (their Fig. 13) normalized to expectations. The latter were, 
in turn, given in terms of ratios to the (tt + + 7r~)/2 yields (AFS Fig. 12). The comparison 
of models to their results would be possible, but would require the recall of information 
from other experiments and the inclusion of the AFS efficiences and acceptance. As we have 
already mentioned, the AFS Collaboration found their data compatible with the other data 
known at the time |L]|| , which will be explored here in detail. We would be unlikely to add 
much to their statement even if we compared || with our models. 

With the data of Antos et al, the situation is different. At first glance they seem to be 
incompatible with all the experiments that have reported the existence of the soft photon 
anomaly. It is therefore very important to include results by Antos et al. in our study. But 
in the paper || they presented the results only in "arbitrary units". A plausible way of 
comparing models to data was suggested to us by J. Schukraft. It is described in the 
relevant subsection below. 
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Unlike in the original paper on the glob model [19 , more detailed and statistically more 



precise results of the simplified version of the glob model can now be combined with the other 
independent sources of photons (MSAM, bremsstrahlung calculated by experimentalists on 
the basis of empirical charged hadron distributions or hadron production models) producing 
a total theoretical expectation. The latter is then compared to the experimental data. We 
can thus also determine the basic parameter of the glob model B = asNcto/Vo more reliably. 



A. BEBC experiment by Chliapnikov et al. [§] 

We start with the BEBC data because they were used in Jl9| to fix the absolute nor- 
malization of the glob model. The experimental cut E l ® b > m n c 2 /2 was enforced in computa- 
tions within both models, as in bremsstrahlung calculations In Fig. ^ we show the data on 
the distribution of photons in x = 2pL,-y/ \fs together with the leading term bremsstrahlung 
calculation from and with the outcome of the glob model. The total theoretical expecta- 
tion, which should be compared to data, is the sum of the glob and bremsstrahlung curves. 
In the glob model calculation we used the value B = 3.5 x 10~ 3 c -1 fm -2 of the overall 
multiplication constant ( |2.5| ). The present value of B, which is 0.7 times the value used in 
|T9]| , was chosen to give good agreement of the theoretical x-spectra not only with the BEBC 



results, but also with the data on prompt photon production in K p and ir + p interactions 
at 250 GeV/c (see below). 

The reader has certainly noticed that the MSAM curve is not shown in Fig. |2|. Neither 
are the MSAM photons included into the total theoretical expectation, which should be 
compared with data. The reason is that their x-distribution is almost flat on the scale 
of Fig. [| and over a wider range it resembles the spectrum of photons from tc° decays 
(not shown). It was the virtue of the experimental procedure used in for isolating the 
anomalous component in the x-spectrum that the photons with a spectrum similar to that 
from the ir° decays were subtracted from the total yield. In this way, the MSAM photons 
were also subtracted and are not present in the data points shown in Fig. ^. 

The transverse momentum spectra of photons produced in 70 GeV/c K + p collisions are 
depicted in Fig. |3|. We can see that the data above px ~ 15 MeV/c are well described 
by the superposition of both models and bremsstrahlung. The glob model provides the 
most important contribution up to pt ~ 45 MeV/c, where MSAM takes over. The classi- 
cal bremsstrahlung formula is dominant below pt ~ 5 MeV/c (as expected from the Low 
formula), but even it is unable to account for all the experimentally observed yield. For 
p T < 10 MeV/c, the mean excess of data above the total expectation (full curve) is roughly 
equal to the total expectation itself. 



B. SLAC experiment by Goshaw et al. 

The inclusive ^-distribution of photons from ir + p interactions at 10.5 GeV/c is shown 
in Fig. The central values of all but one data point lie above the dotted curve, which 
represents a classical bremsstrahlung calculation. But, accounting for experimental errors, 
this does not represent a statistically significant discrepancy. 
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To perform a glob model calculation, we have to rescale the multiplicative constant B to 
the new energy. In [19| it was suggested that it should be done according to the mean hadron 
multiplicity. But we feel that at this very low energy (-y/s = 4.54 GeV) this is inadequate. 
Due to the energy-momentum constraints it is much more difficult to produce a glob with 
mass 1 GeV in addition to a baryon in the final state than a pion (pions account for most 
of the produced multiplicity) with a mass roughly seven times smaller. The dependence of 
Nq on the collision energy should be, at small energies, much steeper than that of overall 
hadron multiplicity and should behave more like, let us say, the mean multiplicity of centrally 
produced 0(1020) 's. Guided by the experimental results on the production of the latter in 
pp collisions (see, e.g., Fig. 3 in [0), we assume that the glob multiplicity at 10.5 GeV/c is 
three times smaller than that at 70 GeV/c. In this way we mimic the threshold effect which 
must exist in the production of globs in low-energy collisions. Scaling according to the ratio 
of produced charged hadron multiplicities in 70 GeV/c K + p and 10.5 GeV/c n + p collisions 
would lead to a decrease by a little smaller factor of 2.1. The ratio of charged multiplicities 
is 1.6. 

In order to be able to compare the model predictions with the data we have taken into 
account the experimental cuts < x < 0.01 and E l ^ b > 30 MeV, number of events (33,676), 
and detection efficiency (0.25) given in Jl|. The results of the glob model are depicted by a 
dashed curve in Fig. || the results of MSAM lie below the lower edge of the diagram. The 
sum of all three theoretical components (solid curve) is now higher than the central values 
of most of the data points, but is as compatible with them as was the pure bremsstrahlung 
component. We can thus conclude that even if the results of the SLAC experiment fl|] do 
not require any additional mechanism rather than classical bremsstrahlung, they do not 
exclude an additional contribution of the size given by the glob model. The glob mechanism 
dominates over bremsstrahlung for px ^ 15 MeV/c. 



C. EMC experiment by Aubert et al. || 

This experiment, which measured the inclusive photon production in deep inelastic scat- 
tering of 200 GeV/c muons in hydrogen, is usually included into a common list with exper- 
iments on anomalous soft photon production in hadronic collisions. We will show here that 
the kinematic range and probably also the production mechanism of the anomaly observed 
in |H are different from those reported by hadronic experiments @J7]||. 

The inclusive photon distributions in || are normalized to all deep inelastic events. We 
will deal with two of them that are presented in transverse momentum px to the virtual 
photon current and fractional energy z, defined as E l ^ b / u, where v is the energy lost by the 
scattered muon. Because we did not have access to empirical distributions in v and in the 
invariant energy of hadronic system W, we fixed them at their mean values (y) = 113 GeV, 
(W 2 ) = 195 GeV 2 , quoted in ||. When calculating the transverse momentum distribution, 
we applied the cut z > 0.05. In virtue of both models we are using here, we assume that 
the intermediate parton systems they are dealing with are bound to the rest frame of the 
hadronic system produced in deep inelastic scattering. 

The ^-distribution of photons from the glob model and MSAM (taken from |TB|) is 



compared with the data in Fig. 61 We can see that the photons from the glob model are 



10 



completely irrelevant here, because their energies are much smaller than those of the prompt 
anomalous photons observed in ||. Having so little energy, the glob photons are below the 
z-cut and do not appear at all in Fig. BL which shows the p\ distribution of anomalous 
prompt photons. The glob model alone provides a successful description of data in both 
cases. 

This points out to a completely different production mechanism of the EMC anomalous 
photons. Rather than being generically related to the anomalous soft photons in hadronic 
collisions, they have more in common with the trimuons discovered a long time ago |36| . It 
is indicated by the fact that also the trimuon production was satisfactorily described by the 



soft annihilation model \Wf 



D. EHS-NA22 experiment by Botterweck et al. @ 

EHS-NA22 Collaboration studied inclusive cross sections of prompt soft photon produc- 
tion in K + p and 7r + p interactions at 250 GeV/c. In our calculations within the glob model 
and MSAM we applied the cut E l ^ b > m 7r c 2 /2, as introduced in the experiment. 

The models are compared with the experimental inclusive photon production cross sec- 
tion in x = 2pi/ ' y/s in Figs. [7| (K + p collisions) and |8] (ir + p collisions). The contribution from 
MSAM is not included for the same reasons as in subsection [IV A| . The agreement between 
data and the combined theoretical expectation is very good. Of course, keeping in mind the 
inevitably large experimental errors (each data point was obtained as a difference of two big 
numbers-the total yield minus the calculated yield from the radiative decays of hadrons), 
theoretical curves scaled down by a factor of, say, ~ 1.5 would be also acceptable. 

In Figs. H and [II], the measured inclusive differential cross sections in transverse mo- 
mentum are presented together with the classical bremsstrahlung estimate made by experi- 
menters, two model curves, and a sum of these three components. For both K + p and n + p 
initial states, the theoretical expectation does not match the data very well. It is below 
the data for ultrasoft photons (jpT ~ 10 MeV/c) and overshoots the data in the medium 
region (10 MeV/c < p T % 50 MeV/c). The former feature is common with other hadronic 
experiments and will be discussed later. The latter may have several origins. First of all, the 



procedure of extrapolations to different collision energies suggested in [|[J| may be unreliable. 
This interpretation is somewhat called into question by the good agreement with the results 
of the WA83/SOPHIE experiment |J (see below), which was done at even a slightly higher 
collision energy {pi a b = 280 GeV/c). Another reason may lie in our inadequate simulation 
of experimental conditions. Besides the energy cut mentioned above, the experimenters 
introduced several others in an effort to minimize systematic errors. The latter would be 
very difficult to implement in our model calculations, because such a task would require a 
detailed knowledge about the experimental setup and a model for the conversion of photons 
to electrons in the metal foils placed inside the Rapid Cycling Bubble Chamber. We are 
not, of course, able to assess here how reliably it was possible to correct for all those cuts, 
and to what extent they persist in the final experimental distributions. 

There is another important issue to be discussed in connection with the paper [[jj. In 
their Figs. 7a,b and 8a,b, the authors of [|7| show the predictions of the glob model, which 
differ a little from what we presented here. In fact, their model histograms are lower than our 
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curves and provide a better description of differential cross sections in p?. Unfortunately, 
the model predictions in [|7| were obtained under oversimplified assumptions. The model 
histograms |19| for K + p collisions at 70 GeV/c were only multiplied by the ratios of total 
inelastic cross sections at 250 and 70 GeV/c. The authors of @ thus neglected the change of 
mean glob multiplicity with the collision energy and assumed that the form of distributions 
is collision-energy independent, x and pt were treated as "scaling variables". 



E. SOPHIE / WA83 experiment by Banerjee et al. [§] 

The high statistics study || of direct soft photon production in n~p collisions at 280 
GeV/c was based on a sample of 310,390 events observed in the apparatus consisting of 
the OMEGA spectrometer and two electromagnetic calorimeters. The results are given 
in the numbers of photons per bin in the variable under consideration (p? or E l ® b ). The 
numbers have been corrected for the 7 detection efficiency, but no attempt has been made to 
correct for the geometrical acceptance and extract the differential cross sections. Generally 
speaking, for experiments with a non-trivial geometrical acceptance that covers only a part 
of the phase space, such a way of presenting results is least biased and therefore most 
valuable Any attempt to go beyond it would require some assumptions about either 
the production mechanism or the photon distribution in inaccessible phase-space regions. 
In order to compare a model with data, one has to modulate theoretical distributions with 
experimental acceptance. 

We took full advantage of the possibilities provided by the simplified version of the 
glob model and merged it with a program that described the geometrical acceptance of the 
WA83 experiment^ in order to calculate the same sort of distributions as shown in ||. Such 
a project would be very difficult to accomplish with a very inefficient photon generator based 
on the original version. For MSAM , the acceptance modulated calculations do not represent 
a novelty pill]. 

The photon transverse momentum spectra from models are compared to the experimental 
one in Fig. |TT]. Again, a reasonable agreement has been achieved for px ^ 10 MeV/c. The 
transition from the glob model to MSAM is now located at pt ~ 30 MeV/c. It is interesting 
that the size of the classical bremsstrahlung is roughly equal to that of the model which just 
dominates. The excess of data over the total theoretical expectation in the ultrasoft region 
(Pt ~ 10 MeV/c) is huge, about sevenfold for the lowest data point. 



F. Soft photon experiment by Antos et al. |J 

As we have already mentioned, the experiment measured the photon production in pBe 
interactions at 450 GeV/c. For those who want to compare its results with theoretical 
estimates, the missing absolute normalization is an obstacle. J. Schukraft, member of both 



I am indebted to M. Spyropoulou-Stassinaki for providing me with the necessary information 
and to A. Belogianni for checking the relevant part of my computer code. 



12 



the HELIOS Collaboration and the experimental group ||, suggested to us to fix the 
model normalization by comparing the bremsstrahlung estimates in || and |J. Because 
of similarities between the two apparatus, they should be identical. But the former is 
given in terms of the double differential cross section, the latter in "arbitrary units". This 
allowed us to recalculate the cross sections provided by models to the "arbitrary units" of 
the experiment 0. For this purpose, we used the Be results of Fig. 13 from the Schukraft 
presentation at the Quark Matter '88 conference M. 

Figure [12] shows the px distribution of photons produced with zero rapidity in the cms of 
proton-nucleon center-of-mass system after subtraction of the decay background. The data 
taken by both detection methods are shown. The systematic errors of the data, the decay 
background and the bremsstrahlung calculation are not reproduced from the original Fig. 
5a to keep our figure uncluttered. The total theoretical expectation agrees nicely with the 
results from the BaF 2 detector and is compatible, taking into account the large systematic 
errors, with the BGO array results. In the kinematic region of the experiment ||, the 
MSAM contributes very little. The glob model provides the most important contibution in 
the medium region 8 MeV/c ^ p-r ~ 40 MeV/c. There is practically no excess in data over 
the theoretical expectation in the ultrasoft region. 



V. CONCLUSIONS AND COMMENTS 

We have seen in the previous section that in the prompt photon production it is useful 
to distinguish among three different kinematical regions: 

(1) ultrasoft, with photon transverse momenta less than, say, 10 MeV/c; 

(2) very soft, 10 MeV/c < p T < 50 MeV/c; and 

(3) soft, characterized by pr ~ 50 MeV/c. 

Of the various components that we included in our phenomenological approach, the classical 
bremsstrahlung formula is invincible in the ultrasoft region. In the very soft region, the main 



contribution to the total theoretical expectation comes from the glob model [19]. The yield 



from the modified soft annihilation model |16| peaks somewhere beyond the range explored 



by the hadronic experiments considered here and is therefore the only theoretical component 
which rises with px- Whether and where it will be above the remaining two components 
(classical bremsstrahlung and glob model) depends on the incident energy and the setup of 
the experiment. But it usually dominates in the soft region. 

The experimental data is mixed, with some experiments showing a much larger excess 
over the hadronic decay background than expected from their classical bremsstrahlung cal- 
culation, and other experiments claiming agreement between the excess and bremsstrahlung. 
The main result of this work is that, for all the experiments the observed x- and pr- distri- 
butions of direct photons are reasonably well described above ~ 10 MeV/c (that is except 
for the ultrasoft transverse momentum region), by a mixture consisting of the classical 
bremsstrahlung calculation (taken from original experimental papers) and two theoretical 
models. The latter are detailed enough to include experimental cuts and acceptances. The 
rules for extrapolating them to different energies have clearly been stated beforehand. On 
the basis of agreement of all data with the theoretical expectation, we conclude that all the 
experimental data (except for the p^-distributions in the ultrasoft region) on anomalous soft 
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photon production are mutually consistent. 

It should be stressed that the magnitudes of different components in the phenomenolog- 
ical mixture did not come out as a result of fitting the experimental data, but are given as 
an interplay between their physics properties and experimental conditions (incident energy, 
instrumental cuts and acceptances). The actual numbers may be very different in different 
cases. For example, in the pt distribution from the BEBC experiment fill , the maximum 



glob/bremsstrahlung ratio is almost seven (see Fig. |^), whereas in the pBe collisions at 
450 GeV/c |J, it does not exceed four (Fig. |12|), and is able to squeeze into the empirical 
upper limits of direct photons provided by the BaF 2 method in || (their Fig. 6). In the 
former experiment, the MSAM gives the dominant contribution for ~ 50 MeV/c, but is 



below the bremsstrahlung up to the highest pxs in Fig. [12]. In the WA83/SOPHIE exper- 



iment ||, the size of bremsstrahlung in the very soft and soft regions is roughly equal to 



that of the dominating model (glob or MSAM, see Fig. |TT 



Let us turn now to the ultrasoft region, where some experiments agree with the theo- 
retical expectation (completely dominated and therefore represented here exclusively by the 
classical bremsstrahlung formula), whereas others see a significant excess. While the deli- 
cacy of the experiments is underlined by the juxtaposition of the early HELIOS || results 
and those of |9[], and the differing results could lie in differing experimental techniques, the 
apparent contradiction between the experiments might also be due to different underlying 
physics conditions. Without any theory or model able to describe the anomalous excess, the 
field is open to speculations. Here is one possibility: 

The data seem to suggest (compare Figs. || || [H| [11], and |l^) that the ultrasoft excess 
decreases with increasing mass of the projectile. There is one experiment that seems to 
contradict this suggestion: the historically first photon experiment [|]], which did not see 
anything anomalous with pions. But let us recall that this is the only experiment which 
used, besides the laboratory energy cut E l ^ b > 30 GeV, also the cut on the longitudinal 
photon momentum in the cms < 2pi iy /y/s < 0.01. The combination of those two cuts 
suppresses the yield of low-p^ photons, especially if they are produced in a narrow cone 
around the projectile momentum. So the absence of the excess in this case need not mean 
its true nonexistence. 

The regularity above may imply that the excess is caused by bremsstrahlung from the 
projectile experiencing a (multiple) small-angle scattering. If we ignore instrumental ef- 
fects (residual gas, thick target, stray photons from background interactions), we can think, 
e.g., about (multiple) soft gluon exchange between projectile and target before the hard, 
multiparticle production interaction takes place. 

Of course, there is still a possibility that the experimenters who saw the excess sig- 
nificantly greater than the classical bremsstrahlung in the ultrasoft region neglected some 
important conventional contribution. 

In the context of the ultrasoft region, it must be also noted that the classical 
bremsstrahlung formula [|l0j, which was used by all experimenters to assess the expected 
level of photon production, is a big unknown. It has been shown a long time ago that it 
should be valid also in the quantum world in the limit of negligibly small photon momenta. 
It enables one to estimate the photon yield if the cross section of the corresponding non- 
radiative reaction is known. For higher photon energies, the non-leading terms in the Low 
expansion p9[ become important. But they cannot be evaluated without a more detailed 
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knowledge of the underlying strong dynamics of the collisions. It is not clear what is the 



region of validity of the classical approximation. In some examples, see, e.g., j40|, it is very 
narrow. 

As already stated in Sec. |, the agreement of the models with experimental data need 
not imply that the mechanisms of photon productions they are based on are real. But let 
us assume for a moment that the models we used in this workP] are more than a clever 
parametrization of all existing data, that they explain the very origin of additional pho- 
tons. We can then go beyond mere phenomenology and address, at least qualitatively, two 
important issues, which may have experimental implications. 

The first remark concerns the dependence of the prompt photon yield on the associated 
hadron multiplicity. The prediction fl2| of faster than linear dependence in dilepton pro- 



duction, and the experiment that seemed to observe it pi] evoked a false impression that 



this effect must take place wherever an anomalous electromagnetic signal is encountered. As 
discussed in more detail in (and, to some extent, already in JT2|), the actual behavior 



depends on the production mechanism. We expect a roughly linear dependence if the dom- 
inant mechanism is bremsstrahlung and faster than linear dependence for the (modified) 
soft annihilation model |12| . In the case of photon production from the glob model, the 



intermediate parton system represents only a part of the event and most of the final state 
hadrons do not originate from it. We therefore expect no correlation between the very soft 
photon production rate and the associated hadron multiplicity. 

Let us also note that a simultaneous observation of two (or more) ultrasoft effects would 
be a nice confirmation of Van Hove's glob mechanism |17|]. It may manifest itself, for 
example, by stronger short-range correlations among pions in events with very soft photons. 

The photon production in high energy collisions remains intriguing and lacking complete 
explanation and therefore deserves continuing experimental and theoretical attention. 
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TABLES 



TABLE I. Parameters of the photon energy distribution in the glob rest frame 



E° (GeV) 


a (fm- 2 GeV- 2 ) 


b (GeV- 1 ) 


c (GeV" 2 ) 


< 0.02 


6.60 x 10 5 


3.39 x 10 1 


1.61 x 10 3 


(0.02, 0.04) 


1.05 x 10 6 


8.15 x 10 1 


4.04 x 10 2 


> 0.04 


1.63 x 10 6 


1.02 x 10 2 


1.82 x 10 2 
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FIGURES 
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FIG. 1. Parametrization of the photon energy spectrum in the glob rest frame (curve) and its 
comparison to the results of original Monte Carlo calculation (crosses). 
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FIG. 2. BEBC data H compared to the classical bremsstrahlung formula (dotted), the glob 
model (dashed), and their sum (solid). 
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FIG. 3. BEBC data [|ll]] and their comparison with the sum (solid) of the classical 
bremsstrahlung formula, the glob model, and the modified soft annihilation model (MS AM). 
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FIG. 4. Goshaw et al. data compared to the classical bremsstrahlung calculation (dotted), 
the glob model (dashed), and their sum (solid). 
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FIG. 5. EMC data pi compared to the modified soft annihilation model (MSAM) and the glob 
model. The rightmost data points can be explained as a bremsstrahlung from the scattered muon 
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FIG. 6. EMC data H compared to the modified soft annihilation model (MSAM). Because of 
the z-cut, the radiation from globs does not contribute at all (compare Fig. ||). 
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FIG. 7. NA22 K+p data compared to the sum (solid) of the classical bremsstrahlung 
culation (dotted) and the glob model (dashed). 
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FIG. 8. Same as Fig. 0, but with a ir + projectile. 
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FIG. 9. NA22 K + p data QJ and their comparison with the sum (solid) of the classical 
bremsstrahlung formula (dotted), the glob model (long dash), and the modified soft annihilation 
model (short dash). 
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FIG. 10. Same as Fig. |, but with a 7r + projectile. 
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FIG. 11. Comparison of the excess in WA83 data || over known sources with the sum (solid 
curve) of the photon yields from classical bremsstrahlung, glob model, and modified soft annihila- 
tion model (MS AM). 




FIG. 12. 450 GeV/c p-Be direct photon data Q compared to the sum (solid curve) of the 
classical bremsstrahlung calculation (straight line), the glob model (long dash), and MS AM (under 
the scale). 
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